Parity Measurement is Sufficient for 
Phase Estimation at the Quantum Cramer-Rao Bound for Path-Symmetric States 
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In this letter, we show that for all the so-called path-symmetric states, the measurement of parity 
of photon number at the output of an optical interferometer achieves maximal phase sensitivity 
at the quantum Cramer-Rao bound. Such optimal phase sensitivity with parity is attained at a 
suitable bias phase, which can be determined a priori. Our scheme is applicable for local phase 
estimation. 
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Interferometry is a vital component of various preci- 
sion measurement, sensing, and imaging techniques. It 
works based on mapping the quantity of interest on to 
the unknown phase of a system and estimating the latter; 
for example, the relative phase between the two modes or 
"arms" of an optical interferometer. Optical interferom- 
etry, often described in the Mach-Zchnder configuration, 
in general differs in the strategies of state preparation and 
detection. The conventional choice is to use a coherent 
light source and intensity difference detection. Assuming 
the unitary phase acquisition operator to be: 



Ua, = e -#(»o-»b)/2 ; 



(1) 



(where h a , fib are the number operators associated with 
the modes,) the phase sensitivity of the conventional 
Mach-Zehnder interferometer (MZI) is bounded by the 
shot noise limit (SNL) 5<f) = \j\ffi (for n photons in the 
coherent state on average); whereas, protocols of quan- 
tum interferometry promise enhanced phase sensitivities 
by prescribing the use of states with nonclassical photon 
correlations and detection strategies that probe the par- 
ticle nature of the output light, for example, via number 
counting As for such nonclassical states, the pro- 
posal to squeeze the vacuum state entering the unused 
port of the conventional MZI, by Caves in 1981, resulted 
in the first instance of a quantum-enhanced MZI capa- 
ble of operating below the SNL |2j. Others, such as the 
twin-Fock state [3j , the maximally path-entangled NOON 
state (\N,Q) + \0,N))/y/2, which reach the Heisenberg 
limit (HL) 6(f) = 1/N, were later proposed 0]. 

Much of the latest experimental efforts in quantum in- 
terferometry have been focussed on attaining the HL [j| . 
The theory of quantum phase estimation aids in identi- 
fying potential schemes for such phase sensitivities Q. 
It is based on the information-theoretical concept of the 
Cramer- Rao bound pj. For the form of unitary phase 
acquisition considered in Eq. ([I]), the quantum Cramer- 
Rao bound (QCRB) — an attribute of the quantum state 
input to the MZI alone (independent of detection) — at 



its best, can reach Heisenberg scaling in the absence of 
photon losses. On the other hand, the classical Cramer- 
Rao bound (CCRB) — an attribute of the combination of 
a quantum state and detection strategy — can optimally 
reach the QCRB of the state. Uys and Meystre de- 
rived optimal quantum states, whose CCRB with number 
counting-based detection attains the HL in the absence 
of photon losses Lee et al. included photon losses 
to the problem and worked out the optimal inputs [9j. 
Meanwhile, Dorner et al. found optimal inputs in the 
presence of photon losses for the generic optimal detec- 
tion strategy based on the symmetric logarithmic deriva- 
tive operation, which is known to achieve the QCRB of 
all quantum states Pezzi and Smerzi revived Caves' 
original scheme, and showed that the QCRB of coherent 
mixed with squeezed- vacuum state input reaches the HL, 
6<f> = l/n, when the two states are mixed in equal pro- 
portions (n/2 photons on average in each state) at the 



input of the MZI 11]. Additionally, they showed that 



the limit could be achieved with the scheme independent 
of the actual value of the phase, by gathering the full 
statistics of number counting in both the modes in place 
of intensity difference detection. 

Apart from number counting and the symmetric loga- 
rithmic derivative operation, another detection strategy 
that has attracted a great deal of attention is the one 
based on the measurement of parity of photon number 
in one of the output modes of the MZI, described by the 
operator II = (—1)™ [12j . Parity detection was first pro- 
posed by Bollinger et al. for enhanced frequency mea- 



surement with an entangled state of trapped ions 13 1. 
Later, Gerry and Campos applied parity measurement 
to optical interferometry with the NOON state for achiev- 
ing phase sensitivities at the HL [LJ]. Parity detection 
achieves sub-shot noise phase sensitivities with various 
inputs [l5[ . MZIs with two- mode squeezed- vacuum state 
input and coherent state mixed with squeezed-vacuum 
state input have been shown to achieve the QCRB of the 
respective inputs with parity detection, which in turn can 



2 



reach the Heisenberg limit [16, 17 1. A theoretical question 
of interest is whether parity detection reaches the QCRB 
of all two-mode input states. In this letter, we prove that 
parity detection achieves the QCRB of all pure states 
that are path- symmetric inside the MZI [18j. The class 
of path-symmetric states includes all pure states that 
have a certain kind of coherence, i.e., phase relation — 
between the probability amplitudes in the photon num- 
ber basis, as given in Eq. (j4}. All single mode state in- 
puts to the MZI result in path-symmetric states after 
passing through the first 50:50 beam splitter. The NOON 
state and the states that result inside the MZI from in- 
puts such as the twin-Fock state [i| , the Yuen state [3] , 
the two mode squeezed-vacuum state [16 L t he coherent 
state mixed with squeezed- vacuum state |ll. 17, 2(3], the 
pair-coherent state 2l|, which are capable of HL phase 
sensitivity, are all path-symmetric states. Based on what 
we prove, it suffices to measure the parity of photon num- 
ber in one of the output modes alone, in place of num- 
ber counting in both the modes, in order to achieve the 
QCRB of all such path-symmetric states, some of which 
can in turn reach the HL. 
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FIG. 1: A MZI with a two-mode input \if)i), which after 
the beam splitter and phase shifter transformations Ubs = 
cxp(-ijJx), U<t> = exp(— i<f)J z ) and U BS = exp(i^Ja;) (in 
that order), is denoted by ^2), 1^3) and \4>a), at the respec- 
tive stages. 



We work in the Schwinger representation, wherein two- 
mode interferometry is described in terms of SU(2) al- 
gebra [22j. A two mode A-photon state in this repre- 
sentation resides in the j — N/2 subspace of the angu- 
lar momentum Hilbert space. A typical MZI, Umzi = 
U bs U^UbSi as described in Fig. [TJ is considered. For 
the chosen phase shifter and beam splitter transfor- 
mations Utf, — exp(—i(/>J z ) and Ubs — exp(— i-fJ x ), 
Umzi = exp(— i</>J y ). We call this the type-I MZI. On 
the other hand, if we choose a beam splitter that per- 
forms Ubs — exp(— i\ J y ), then Umzi = cxp(— i(f>J x ), 
which we will call a type-II MZI. Without loss of gen- 
erality, one can switch from type-I to type-II MZI by 
suitably adding a finite number of phase shifters before 
and after the beam splitters. 

The error in the estimate of the unknown phase shift 
d> based on the measurement of an observable O can be 



written based on the linear error propagation formula as 

5(f) = AO/(d(0)/d(p) . For an observable that acts on 

the state ^3) , this quantity S(f>, due to the Heisenberg un- 
certainty principle between the generator of phase evolu- 
tion J z and the observable, obeys (8(f) > 1/(2AJ 2 ). For 
pure quantum states \ip3), the right-hand side of this in- 
equality is a tight bound, and is identically equal to their 
QCRB. The equivalent (necessary and sufficient) condi- 
tion for achieving the bound is that of equality in the 
Heisenberg uncertainty relation, namely: 



[0-{0)I)\^) = i\{J z -{J z )I)\^), 



(2) 



where A is purely real and the expectation values are 
calculated with respect to the state |V>3)- 

We first examine photon number counting-based detec- 
tion strategies. For mathematical convenience, we choose 
a type-II MZI. Since photon counting at the output of a 
MZI is equivalent to the measurement of the operator J z 
on the state IV^}; the corresponding observable O acting 
on ^3) for such detection strategies must be diagonal 
in the J x basis. This is so, because the J z basis states 
transform into the J x basis states under the action of the 
type-II beam splitter transformation Ubs — ex p(~ i^Jy)- 
In the J x basis {Ito^)}, since the matrix elements of J z 
are purely imaginary, the right hand side of Eq. @ is 
real (Hermiticity requirement for the observable O) if 
and only if the coefficients of 1 1^3) are purely real (up 
to a global phase factor), i.e.: 

(m x \^) = (m B |V3)*e- Q x, V m x € (3) 

where x is a real constant, and ( J z ) = 0. This condition, 
when written in the J z basis {|ro z )}, transforms into: 

(mM = (-m^ 3 )*e-* 2 * V m z e {-j, +3}. (4) 

(Note that ( J z ) of such states is always zero.) There- 
fore, all pure states that obey Eq. ((4]) reach their maxi- 
mal phase sensitivities with detection strategies based on 
photon number-counting. (For the rest of the paper, we 
will use til = m z .) 

Further, if |V> 2 ) = J2m=-j C ™I TO ) > then IV'3) = 
Sm=-7 c m e~ lmct '\m) , which implies: 



(ra|V>3 



(m\i> 2 



(-m|V> 3 >* (c_ m e^™)* c*_ m {-m\i) 2 y 



(5) 



However, (m\ip^) / (— m\ip3)* = e~ t2x . Therefore 
(m\tp2) / {—m\tp2}* also equals e~ l2x , i.e., the condition 
of Eq. dU can be satisfied independently of the value of 
the unknown phase (p. In a recent paper 18], Hofmann 
interpreted this phase independence of Eq. (j4} as a sym- 
metry in the Heisenberg picture, wherein, the angular 
momentum vector J = {J x , J y , J z } associated with the 
state is invariant under rotations about the z-axis 231. 
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He called such states as path- symmetric states, since in 
the Schrodinger picture, the condition in Eq. ((4]) can be 
interpreted as the invariance of the two-mode pure state 
\ip ) with respect to a mode exchange operation (along 
with complex conjugation of the coefficients). 

We now study such path-symmetric states in the con- 
text of parity-based metrology. Parity has been com- 
monly studied with MZIs containing the beam split- 
ter transformation Ubs — exp(— i\J x )\ therefore we 
switch to the type-I MZI. The observable O (at the state 
\ip 3 )) corresponding to parity measurement in the output 
mode bf can be obtained by transforming II = ( — l)i~ Jz 
through the beam splitter. It is given by the operator 
Q = UbsWIs M- 



N 

aN \ " 



(-l) k \k,N-k){N-k,k\ 



k=0 
3 



J2 (-iy +m n(-m\, 



(6) 



(7) 



and (-04|n|-04) = (ip 3 \Q\ip 3 ). In order to achieve the 
quantum Cramer-Rao bound with the operator Q for a 
path-symmetric state \ip 3 ), we need: 



(8) 



where (Q) = ^ ^ Ar (-l) J+m c_ m c I *„e ^2m * and A = 
m=—j 

±[AQ/AJ Z ] is a real number. Since Q 2 = I, alterna- 
tively, we have 



(m\Q\ip 3 ) - {Q)(m\i/) 3 ) 



±i 



1 - (Q? 



{J: 

1 - (Q? 



1/2 



{m\J z \ip 3 ) 

Via direct computation of each term, we obtain: 

1/2 

I — ( ): - 

S — (Q) = ±im 



(9) 



(10) 



(J* 

where S = ^(_i)J+™^e i2m * = S' + iS". Let c r , 
T m e %6m , where r m and 9 m are real. Then, 

5" = Re (i N '(_i)J+™ e o("»*-x-0m) 

_ j±(-l) j+m sm2(m(t)-x-0 m ) if N is odd, 
~ \{~\) m cos2{mct)~x-0m) if N is even. 

(11) 

If S' = ±1 for some cf> and for all to, then S" = 

2 



since S' z + S' u = 



that case, Re (i N {-l) J+m c- m c* m e 



i N (— i)J+ m e i2 ( m ^ _ x-0m) 



* i2m<fi\ 



= 1. In 



Im (^(-1)^+' 

(Q) = E ^(-1) 

m=-j 



V r* P i2rn 4 



= for all m. So 



j-\-m 



* Alrad 



and 



±1, 



which implies that the condition in Eq. (|T0f is satisfied, 
and hence the quantum Cramer-Rao bound is achieved. 
Therefore, S' = ±1 is a sufficient condition for a path- 
symmetric state to achieve its quantum Cramer-Rao 
bound with parity measurement. 

Let us assume that phase sensitivity at the quan- 
tum Cramer-Rao bound is obtained for a path-symmetric 
state IV'2); i-e., the equation (|10[) holds for all to; 



S' + iS" - (Q) = ±im 
Then, we have 

S' = (Q), S" = ±m 
and furthermore, 



i-(Q) ; 



i-(Q) ; 



1/2 



(12) 



1/2 



(13) 



5'" 



±to 



1 



(Jz 



1/2 



±mS" 



(Jz 



1/2 



(14) 



If \S'\ ^ 1, then S" ^ since S' 2 + S" 2 = 1. So we get 

(Jz ) ly/2 = ±ti for all to, but this result is a contradiction 
in general, since it should hold for all to. Hence, |S"| = 1. 
Therefore, S' = ±1 is also a necessary condition for a 
path-symmetric state to achieve its quantum Cramer- 
Rao bound with parity measurement. Hence, we say 
an N -photon path- symmetric state \1fi2} = zCm=-j c m|™) 
reaches its quantum Cramer-Rao bound with the parity 
detection operator Q if, and only if, there exists a real 
number <j) such that, for all m, 



Re (V(-l) J ' +m ^e i2m < 



±1. 



(15) 



Although seemingly complicated at the outset, the con- 
dition in Eq. (p~5)) is easily satisfied in the following cases: 

(a) c m = r m , i.e., c m is purely real (9 m — 0) for all to, 
and therefore x = 0. 

(b) c m = r m e" ne , i.e., 9 m = m9 for some constant 9, 
and therefore \ = 0. 

The values of <f> in case (a), and <f) — 6 in case (b), when 
chosen as: 



if N is odd, 



— nr 3^ 
2 U1 2 

or ^ or 7T or ^ if is even, 



m=-j 



(16) 



satisfy the condition in Eq. (115]) , and consequently re- 
sult in maximally sensitive phase estimation with parity 
measurement for the given state at those values of </>. 

Additionally, given any arbitrary path-symmetric state 
^2) with real parameters {r m ,9 m }, and a real constant 
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X, there always exists a bias phase (3, such that 9 m +mf3 = 
— x for all m; and the value of /3 can be obtained from: 



Y J {e m + m(3) = -Y j 



(17) 



m>0 



m>0 



Solving Eq. (ITTl) . we get: 



/3 = 



_8 ^(JV+l) 
(N+l) 2 V 2 A 



E 



m>0 







' AT(JV+2) \ 2 



m>0 







if N is odd, 
if N is even. 

(18) 

in Eq. (fTTj) with m/3, and choos- 
as per Eq. (HU), the QCRB of the 



By replacing — x - 
ing the value of <fi 
input state can be reached with the operator Q at the 
corresponding value of (f>. So, for any given arbitrary N- 
photon path symmetric state, there always exists a real 
phase 4>, wherein parity detection attains the QCRB of 
the state. These results directly extend to states with 
indefinite photon numbers, with the optimal point of op- 
eration (the value of <fi or <fi — 9 or <fi — (3, as is the case) 
given by the values common to the even and odd compo- 
nents of Eq. (16]), i.e. (21 + 1)tt/2, I e Z. 

Although optimal phase sensitivities with parity are 
reached only at specific bias phases, yet, since the bias 
can be predetermined, parity detection provides a use- 
ful tool for "local" phase estimation. The biases may be 
easily implemented in the laboratory if we notice that 
the phase <fi is the relative phase between the two arms 
of the interferometer, and hence can be always tuned to 
the optimal bias point or "sweet spot", by inserting a 
an adjustable and known phase shift in the other arm. 
For detecting small changes of more or less known pa- 
rameters over space or time, a priori knowledge about 
the parameter is assumed — this is called local parameter 
estimation — as opposed to the "global" one, wherein a 
complete lack of knowledge about the parameter is ini- 
tially assumed [24j. 

One way to measure parity is obviously to perform 
number counting at the output and to infer the par- 
ity from it. High precision number-resolving detection 
of photons at single-photon level has been demonstrated 
experimentally with the superconducting transition-edge 
sensors [25j |. Apart from number counting, there have 
been other methods proposed for the direct measurement 
of parity. Gerry and co-workers suggested the use of opti- 
cal nonlinearities 14J . Plick et al. showed that homodyne 
quantum state tomography can be used to construct the 
expected parity signal, at least in the case of Gaussian 
states, since the expectation value of the parity operator 
is proportional to the value of the Wigner function of the 
state at the origin in phase space for such states [26| . 

In conclusion, we have shown that the detection strat- 
egy based on the measurement of the parity of photon 
number achieves the maximal phase sensitivity at the 
quantum Cramer-Rao bound of all pure states that are 



path-symmetric. Although such optimal phase sensitivi- 
ties with parity detection are achieved only locally, those 
specific bias phases can be predetermined. 
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tional Science Foundation. 
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